Improved Thermoelectric Cooling Based on the Thomson Effect 
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Traditional thermoelectric Peltier coolers exhibit a cooling limit which is primarily determined by 
the figure of merit, zT. Rather than a fundamental thermodynamic limit, this bound can be traced 
to the difficulty of maintaining thermoelectric compatibility. Self-compatibility locally maximizes 
the cooler's coefficient of performance for a given zT and can be achieved by adjusting the relative 
ratio of the thermoelectric transport properties that make up zT. In this study, we investigate the 
theoretical performance of thermoelectric coolers that maintain self-compatibility across the device. 
We find such a device behaves very differently from a Peltier cooler, and term self-compatible 
coolers "Thomson coolers" when the Fourier heat divergence is dominated by the Thomson, as 
opposed to the Joule, term. A Thomson cooler requires an exponentially rising Seebeck coefficient 
with increasing temperature, while traditional Peltier coolers, such as those used commercially, have 
comparatively minimal change in Seebeck coefficient with temperature. When reasonable material 
property bounds are placed on the thermoelectric leg, the Thomson cooler is predicted to achieve 
approximately twice the maximum temperature drop of a traditional Peltier cooler with equivalent 
figure of merit (zT). We anticipate the development of Thomson coolers will ultimately lead to solid 
state cooling to cryogenic temperatures. 

PACS numbers: 84.60. Rb, 05.70.Ce, 72. 20. Pa, 85.80.Fi 



1. INTRODUCTION 

Peltier coolers are the most widely used solid state 
cooling devices, enabling a wide range of applications 
from thermal management of optoelectronics and infra- 
red detector arrays to polymerase chain reaction (PCR) 
instruments. Thermoelectric coolers have been tradition- 
ally understood by means of the Peltier effect, which de- 
scribes the reversible heat transported by an electric cur- 
rent. This effect is traditionally understood in terms of 
absorption or release of heat at the junction of two dis- 
similar materials. The conventional analysis of a Peltier 
cooler approximates the material properties as indepen- 
dent of temperature (Constant Property Model (CPM)). 
This results in a maximum cooling temperature differ- 
ence AT max for a CPM cooler, which dependent on the 
figure of merit ZT of the device [1, 2]. 

AT max = (1) 

For the best commercial materials this leads to a 
AT max of 65K (single stage) [3], which translates to a 
device ZT at 300K of 0.74. In the CPM the device ZT 
is equal to the material zT. Material zT depends on the 
Seebeck coefficient (a), temperature (T), electrical resis- 
tivity (p), and thermal conductivity (n), zT — a -^-. In 
the CPM, the only way to increase AT max for a single 
stage is to increase zT, leading to the focus of much ther- 
moelectric research on improving zT. It is well known 
that even further cooling to lower temperatures can be 
achieved using multi-stage Peltier coolers [1, 2]. In prin- 
ciple, each stage can produce additional cooling to lower 
temperatures, regardless of the zT of the thermoelectric 



material in the stage. In practice, the thermal losses 
and complications of fabrication limit the performance 
of such devices. The 6-stage cooler of Marlow achieves 
a AT max of 133 K; this doubling of AT max compared to 
a single stage cooler is achieved despite using materials 
with similar zT [3]. Alternatively, such AT max with a 
single-stage CPM cooler would require ZT to be 2.5. 

The transport properties across a single thermoelec- 
tric leg can be manipulated to improve cooling perfor- 
mance, although it has been less effective in reducing 
AT max than a multi-stage approach. One common strat- 
egy is to engineer a change in extrinsic dopant concen- 
tration across a thermoelectric element which can sig- 
nificantly alter a, p and even k. For example, this has 
been demonstrated for thermoelectric generators in n- 
type PbTe doped with I [4]. Similar efforts have been 
done with cooling materials, as has been reviewed in ref 
[ ]. The simplest explanation for an improvement is an 
in increase in the local zT at some temperatures by spa- 
tially adjusting the dopant composition within a material 
[6]. 

Early theoretical work by Sherman et al for TEC found 
that different AT max could be predicted from materials 
have the same or similar average zT but different temper- 
ature dependence of the individual properties a,p,K [ ]. 
This demonstrated that optimizing cooler performance is 
significantly more complex than simply maximizing zT. 
More recently, Muller et al. [8-10] and Bian et al. [11, 12] 
used different numerical approaches to predict substan- 
tial gains in cooling to AT max from functionally grading 
where an average zT remains constant in an effort to 
determine the best approach to functionally grading. 

Different material classes optimized for different tern- 
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Fig. 1: The local reduced coefficient of performance <j> r is 
optimized at a specific reduced current density, termed s. If 
u/s, the (j> r is less than that predicted by the material zT. 
Here, z = 0.002 K~\ a = 200/i^A'- 1 . 



peratures can also be segmented together to improve per- 
formance of thermoelectric generators but the current 
must also be matched [13]. The analysis of segmenta- 
tion strikingly demonstrates that increasing the average 
zT does not always lead to an increase in overall thermo- 
electric efficiency and so an understanding of the thermo- 
electric compatibility factor is needed to explain device 
performance [ ]. 

This paper derives the cooling limit for a single stage, 
fully optimized (self-compatible) TEC that functions as 
an infinitely staged cooler. The Fourier heat divergence 
in such an optimized cooler is found to be dominated 
by the Thomson effect rather than the Joule heating 
as in traditional Peltier coolers. This new opportunity 
presents a new challenge for material optimization based 
on compatibility factor rather than only zT . 



2. THEORY 

Coolers are characterized by the coefficient of perfor- 
mance (cj) = Qc/P), which relates the rate of heat extrac- 
tion at the cold end Q c to the power consumption P in 
the device [15]. For simplicity, but without loss of gener- 
ality, a single thermoelectric element can be considered 
rather than a complete device. A TEC leg can be treated 
as an infinite series of infinitesimal coolers, each of which 
is operating locally with some COP. Scaling this COP to 
the local Carnot COP (T/dT) yields the local reduced 
coefficient of performance <j) r . [16]. This relationship be- 
tween local performance across the leg and global COP, 
(j>, given in Eq. 2 is derived in the Appendix based on 
Ref.[7][17]. While TECs are traditionally analysed using 
a global approach, we have previously shown the utility of 
a local approach [14, 15, 18]. This local approach leads to 
a consideration of material 'compatibility', as discussed 
below. 
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Fig. 2: a) The CPM Peltier cooler and u — s Thomson cooler 
are compared using the same constant z = 0.002 K _1 . The 
overall device <j> of a CPM cooler crosses zero at a finite tem- 
perature, indicating AT max is reached, while <f> remains pos- 
itive for all temperatures for the u = s cooler, b) The local 
performance of a CPM cooler (cj> r ) is significantly compro- 
mised at both the hot and cold ends. In contrast, 4> r ,max is 
achieved at all temperatures when u = s. In both panels, the 
performance calculated for an actual Bi2Te3 Peltier cooler leg 
is similar to the CPM. 

The compatibility approach to optimizing thermoelec- 
tric cooling arises naturally from an analysis of the ther- 
mal and electric transport equations. This method has 
been described in detail for thermoelectric generators [15] 
and coolers[19] and are reproduced here for TEC. The 
method has been experimentally verified [2( ] and shown 
to reproduce results using a more traditional finite el- 
ement results but with less computational complexity. 
This method has been incorporated into several engineer- 
ing models such as those used by NASA for Radioisotope 
Thermoelectric Generators [21] and Amerigon/BSST for 
automotive applications [20, 22]. Consider an infinites- 
imal section of thermoelectric leg in a temperature gra- 
dient and an electric field. The temperature gradient 
will induce a Fourier heat flux (q K = — kVT) across this 
segment. The divergence of this heat (Eq. 3) is equal to 
the source terms: irreversible Joule heating {pj 2 ) and 
the reversible Thomson heat (Tj^jVT), both of which 
depend on the electric current density (J). From these 
two effects, the governing equation for heat flow in vector 
notation is 

V • q K = V ■ (-«VT) - pj 2 - t j • VT (3) 

with Joule heat per volume pj 2 , Thomson coefficient 
t = and Thomson heat per volume r j • VT. The 
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Peltier, Seebeck and Thomson effect are all manifesta- 
tions of the same thermoelectric property characterized 
by a. The Thomson coefficient (r = Tj^) describes the 
Thomson heat absorbed or released when current flows 
in the direction of a temperature gradient. 

Restricting the problem to one spatial dimension, Eq. 
(3) is typically examined assuming the heat flux and elec- 
tric current are parallel [15]. In the typical CPM model 
used to analyze Peltier coolers, the Thomson effect is zero 
because a is constant along the leg =0). 

The exact performance of a thermoelectric leg with 
a(T), p{T), and k(T) possessing arbitrary temperature 
dependence can be straightforwardly computed using the 
reduced variables: relative current density (u) and ther- 
moelectric potential ($) [ ■]. The relative current den- 
sity u, given in Eq. 4, is primarily determined by the 
electrical current density j, which is adjusted to achieve 
maximum global COP. The thermoelectric potential $ is 
a state function which simplifies Eq. 2 to Eq. 6 [15]. 
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(5) 



(6) 



Changing variables to T via the monotonic function 
x(T), Eq. 3 simplifies to the differential equation in u(T). 



du 
dT 



u 2 I T 



da 
dT 



(7) 



The maximum local 4> r , denoted 4> r ,max, occurs when 
u = s. The expression for 4>r,max (Eq. 10) is an explicit 
function of the material zT and is independent of the 
individual properties a, p, n. This maximum allowable 
local efficiency provides a natural justification for the def- 
inition of zT as the material's figure of merit. 



Vl + zT - 1 
Vl + zT+l 



(10) 



One thus wishes to construct devices where, locally, each 



segment has "it = s" and thus 



is obtained. Glob- 



ally, maximum <j) is found when the entire cooler satisfies 
u = s. 



3. COOLING PERFORMANCE 

To compare the cooling performance of traditional 
Peltier coolers and u = s coolers, we consider coolers 
with equivalent z. Traditional Peltier coolers have typ- 
ically been analyzed with the constant property model 
(CPM), yielding a constant z (where zT is linearly in- 
creasing with temperature). We will show that constant 
z, but allowing a, K, p to vary with T, can lead to sub- 
stantial improvement in cooling. At the limit of this 
variation, we will assume the properties can be varied 
to satisfy u = s. 

Performance of a CPM cooler CPM coolers have 
been extensively studied, typically using a global ap- 
proach to the transport behavior. The <fr for a CPM 
cooler (operated at optimum j) is given by Eq. 11 [24]. 
Figure 2a shows the of a CPM cooler decreases with 
increasing AT. With increasing cooling, this cj> decreases 
and reaches zero at AT max (Eq. 1). 



Using this formalism, the reduced coefficient of perfor- 
mance (<p r ) can be simply defined for any point in the 
cooler (Eq. 8). Fig. 1 shows this relationship between u 
and 4> r . From Eq. 2, it can be shown that <j> is largest 
when <j) r is maximized for every infinitesimal segment 
along the cooler. Hence, global maximization can be 
traced back to local optimization [23]. 
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(8) 



The optimum u which maximizes cf> r {-P 1 - = 0) can 
be expressed solely in terms of local material properties 
(Eq. 9). This optimum value of u is defined as the ther- 
moelectric compatibility factor s c for coolers. 



VI + zT-l 
aT 



(9) 



As this paper strictly focuses on coolers, we will refer to 
s c as simply s. 



lCPM 




(11) 



To understand what is limiting the CPM cooler at 
&T max , we derive the local reduced coefficient of per- 
formance <j)^ PM {T). To obtain (j)^ PM we need u as a 
function of T. The solution to differential equation 7 for 
CPM is 
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= -f + (T h -T) 



(12) 



where the value of u at T = Th (uh) serves as an ini- 
tial condition. This expression allows u(T) to be de- 
termined for any CPM cooler, regardless of temperature 
drop (AT < AT max ) and applied current density (j). 
The global maximum COP (4>) is obtained when the op- 
timum Uh from Eq. 13 is employed. 
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Consideration of Eq. 13 reveals that the maximum T c 
is obtained when 1/tt/j approaches zero. Figure 3 shows 
|u| becoming infinite at Th for the CPM cooler. In this 
limit, Eq. 13 can be simplified to give Eq. 1 with Z = z. 
Thus, a local approach to transport yields the classic 
CPM limit typically obtained through an evaluation of 
global transport behavior. 

Combining Eq. 8, 12, and 13 results in 4> r (T) at AT max 
for the CPM Peltier cooler (Eq. 14). This expression re- 
veals 4>r drops to zero at both ends of the CPM cooler leg, 
as shown in Figure 2b. This prohibits additional cooling 
and sets AT max . 
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To achieve cryogenic cooling (T c — > 0) within the CPM, 
zT must approach infinity (Eq. 1). For example, cooling 
with a single-stage CPM cooler to 10 K would require 
zT to be over 1000 if the hot side is 300 K. When is 
negative, the net effect of the thermoelectric device is to 
supply heat, rather than remove heat, from the cold side. 
For negative (j) values for the CPM cooler to be obtained 
requires certain parts of the cooler to locally possess (j> r < 
0. Such a result may be surprising at first as this (f> r < 
region is made from material possessing positive zT. This 
seems particularly odd when compared to the behavior of 
staged generators, discussed above. Clearly, single- and 
multi-staged CPM legs exhibit fundamentally different 
behavior, despite being composed of exactly the same 
material. Such behavior can be rationalized using the 
thermoelectric compatibility concept. 

Figure 3 shows that the compatibility condition (u = s) 
is maintained at only one point in the CPM cooler. Con- 
sequently, CPM coolers operate inefficiently (u 7^ s) at 
both the hot and cold ends. This is demonstrated in Fig- 
ure 2b, where <fi r < 4>r,max for ah but one point. Once 4> r 
goes below zero at low temperature, the thermoelectric 
device is no longer cooling the cold end and AT max is 
reached (Figure 2a). 

While real coolers do not possess temperature- 
independent properties, the qualitative results for CPM 
translate well to traditional Peltier coolers due to their 
weak material gradients. Considering a Bi2Te3 leg with 
temperature-dependent properties described in Ref. [25] , 
we find u and s to be quite close to a z-matched CPM 
cooler (Figure 3). Like the CPM cooler, u = s at only one 
temperature along the leg. This leads to similar (j) r {T) 
for the Bi 2 Te3 and CPM coolers, shown in Figure 2b. 
Within the CPM, large zT results in a high upper 

, ax is achieved. 



limit to 



but does not ensure this 



Generaly, commercial cooling materials such as Bi2Te3 
and any material that can be described by the CPM 
model will be operating significantly below the <fi r 
predicted by the zT they possess (Eq. 8, 10). 

Performance of a u — s cooler We now consider an 
idealized cooler which maintains u — s across the entire 



Fig. 3: Both the CPM and Bi2Te3 coolers have u = s at one 
point along the leg. By definition, the u = s is self-compatibile 
along the entire leg. The different slope signs of u for CPM 
and u — s reveals that these coolers are fundamentally dis- 
tinct. The curves were generated for an optimized cooler at 
ATmax with z = 0.002 KT\ T h = 300 K. 



leg. 4> r for this cooler is simply given by Eq. 10. This <p r 
is found to be positive for all T, as z is always a posi- 
tive real number. Globally, this translates to the analytic 
maximum for <f> for a cooler where z is defined and lim- 
ited. 

To facilitate comparison with CPM, we consider a con- 
stant z model where the individual properties are ad- 
justed to maintain u = s. The constant z approach 
yields vanishing zT at low T, consistent with real ma- 
terials. Evaluating tfi (Eq. 2) for a u = s cooler and 
the assumption of constant z, one obtains Eq. 15, where 
Mi = y/T+zTi with Ti = T h , T c . 



(M h 



exp 



2(M h - M c ) 
(M h - 1)(M C - 1) 



-1 (15) 



Inspection of Eq. 15, where Mh > M c > 1, reveals that <fi 
is always greater than zero for a u = s cooler. 

The difference between CPM and u — s coolers can 
be visualized in Fig. 2a, with the of the Thomson 
cooler asymptotically approaching zero with increasing 
AT. Figure 2b shows that <j> r for a self-compatible cooler 
with constant z remains finite and positive throughout 
the device. In contrast, the CPM cooler is operating 
inefficiently at both the hot and cold ends, limiting its 
temperature range. 

In principle, if u = s can be maintained, the idealized 
u = s cooler can achieve an arbitrarily low cold side 
temperature as long as the all of the materials have 
a finite zT. However, the material requirements to 
maintain u = s become exceedingly difficult to achieve 
as the cooling temperature is reduced and the ultimate 
cooling will be finite, yielding T c > 0. 
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4. MATERIAL REQUIREMENTS 

A CPM cooler has a fixed z and performance which 
is independent of the ratio of individual properties as 
long as they are constant with respect to temperature 
(Eq. 1). In contrast, a u = s cooler requires dramatic 
changes in properties with temperature to maintain self- 
compatibility. Within the constraint of constant z, con- 
sideration of Eq. 9 suggests that the Seebeck coefficient 
must be varied across the device to maintain u — s. Ad- 
ditionally, as a (T) = y/zp(T)K(T) within a constant z 
model, the product p(T)n(T) must also vary across the 
device. 

The Seebeck coefficient profile a(T) for a u = s cooler 
with constant z can be solved analytically. Combining 
Eq. 7 and u = s yields the simple differential equation of 
a(T): 



d ( aT \ _ da a 1 + %/l + zT 



Solving this equation yields 



a(T)= ao ^_^ex P 



Vl + zT- 1 



(17) 



With this expression for a(T), it is possible to eval- 
uate s(T) with Eq. 9. Figure 3 shows the variation in s 
required for a u — s cooler with constant z. The self- 
compatible cooler modeled in Figure 3 has z — 0.002; 
60 K of cooling results in a change in s of one order of 
magnitude. 

The approximation for small zT yield a simple expres- 
sion for a(T), given by Eq. 18. 



dT 



(lna(T)) 



4 

Tr 2 



a(T) oc exp 



-4 

zf 



(18) 



This reveals that a should be very large at the hot end 
and must decrease to a low value at the cold end. This 
exponentially varying a(T) required to maintain u = s 
for constant z is anticipated to be the limiting factor in 
real coolers and place bounds on the maximum cooling 
obtainable. We consider the realistic range of a below. 

Large values of a are found in lightly doped semicon- 
ductors and insulators with large band gaps (E g ) that 
effectively have only one carrier type, thereby preventing 
compensated thermopower from two oppositely charged 
conducting species. Using the relationship between peak 
a and E g of Goldsmid (Eq. 19) allows an estimate for the 
highest a(Th) we might expect at the hot end, ah [26]. 
Good thermoelectric materials with band gap of 1 eV are 
common while 3 eV should be feasible. For a cooler with 
an ambient hot side temperature, this would suggest ah 
should be ~l-5mV/K. Maintaining zT at such large a 
will require materials with both extremely high electronic 
mobility and low lattice thermal conductivity. 



A lower bound to a c also arises from the interconnected 
nature of the transport properties. We require zT to be 
finite; thus the electrical conductivity a must be large as 
a c tends to zero. In this limit, the electronic component 
of the thermal conductivity {ke) is much larger than the 
lattice (k_l) contribution and k ~ k^. To satisfy the 

2 > 2 

Wiedemann- Franz law (ke = LaT where L — is 
the Lorenz factor in the free electron limit), a c has a 
lower bound given by Eq. 20. For example, a z = ^ 
K _1 and T c = 175 K results in a lower bound to a c of 
119/xV/K. 



LzT r 



-zT r 



(20) 



The maximum cooling temperature T c can be solved as 
a function of z, E g and Th from equations Eq. 17, Eq. 19 
and Eq. 20. For small z the approximate solution 
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(21) 



gives an indication of the important parameters but 
quickly becomes inaccurate for zT above 0.1. 

Material limits to performance With these bounds on 
material properties, we consider the AT. max of a it = s 
cooler. Figure 2 suggests that the of a u = s cooler 
remains positive for all temperature. However, obtaining 
materials with the required properties limits AT max to a 
finite value. Fig. 4 compares the AT max solution for u — 
s and CPM coolers with the same z. Here, the maximum 
Seebeck coefficient is set by the band gap (E g = 1 — 3 eV), 
per Eq. 19. The u = s cooler provides significantly higher 
AT max than the CPM cooler with the same zT, nearly 
twice the AT max for E g — 3 eV. 

Spatial dependence of material properties These an- 
alytic results are possible because the compatibility ap- 
proach does not require an exact knowledge of the spatial 
profile for the material properties. Nevertheless, it is pos- 
sible solve for the spatial dependence of the u = s cooler, 
given some material constraints. To determine x(T), we 
integrate Eq. 4, recalling we have assumed constant cross- 
sectional area (j(x) = const.), obtaining Eq. 22. 



-1 [ ±h 
x(T) = —r- / undT 
3 Jt 



(22) 



a h = E g /{2eT h ) 



(19) 



Thus, the natural approach to cooler design within the 
u = s approach is to determine the temperature depen- 
dence of the material properties, and then determine the 
required spatial dependence from the resulting u(T) and 
«(T). 

Figure 5a shows an example of the Seebeck distribution 
a(x) along the leg that will provide the necessary a(T), 
where a constant kl =0.5 W/mK is assumed. The a of 
Figure 5a spans the range permitted by Eq. 19 and 20. 

In a real device the spatial profile of thermoelectric 
properties will need to be carefully engineered. If this 
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rapidly changing a(x) is achieved by segmenting differ- 
ent materials, low electrical contact resistance is required 
between the interfaces. We anticipate such control of 
semiconductor materials may require thin film methods 
on active bulk thermoelectric substrates. 
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Fig. 4: The maximum temperature drop AT m ax of a u — s 
Thomson cooler exceeds that of a Peltier cooler with the same 
- Large band gap, E g , thermoelectric materials are necessary 
the hot junction improves the performance (2), = 300K). 
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Fig. 5: a) The Seebeck coefficient of the u — s Thomson 
cooler varies exponentially, while it is by definition constant 
for the CPM cooler. As a degenerately doped semiconductor, 
the Seebeck coefficient of commercial Bi2Te3 increases gradu- 
ally with temperature, b) The curvature of T(x) for the CPM 
Peltier cooler temperature profile is opposite that of the u = s 
cooler because of the different sign of the Fourier heat diver- 
gence. Again, similar behavior is found between the CPM 
and commercial Bi2Te3 cooler. 



5. COOLER PHASE SPACE 

The improved performance of a u — s cooler is not 
simply an incremental improvement, but rather we find 
CPM and u — s coolers operate in fundamentally differ- 
ent phase-spaces. Here, by phase space we refer to the 
class of solutions defined by the sign of the Fourier heat 
divergence (V ■ q K in Equation 3). The Fourier heat di- 
vergence in a cooler contains both the Joule (pj 2 ) and 
Thomson (rj • VT) terms. 

We begin by considering the Fourier heat divergence in 
CPM and Bi 2 Te3 coolers and then compare this behavior 
to ii = s coolers. In the typical CPM model to analyze 
Peltier coolers, r = as there is no variation in a. In a 
CPM cooler, V • q K is thus greater than zero. This can 
be seen by the downward concavity of the temperature 
distribution in Figure 5b. In a typical Peltier cooler (e.g. 
Bi2Te3), the concavity is the same as the CPM cooler and 
thus the divergence is likewise positive. This is because 
the Thomson term is always less than the Joule term in 
a conventional thermoelectric cooler. 

In contrast, a u — s cooler changes the sign of the 
Fourier heat divergence such that V • q K is less than zero. 
This can be readily visualized in Figure 5b, where the 
concavity of the u = s cooler is opposite the CPM and 
Bi2Te3 coolers. This difference in concavity must come 
from the Thomson term being positive and greater than 
the Joule heating term. The large magnitude of Thom- 
son term is understandable with the exponentially rising 
Seebeck coefficient seen in Figure 5a. The reversibility of 
the Thomson effect requires that for V • q K to be less than 
zero, the hot end must have a high \a\ relative to the cold 
end, and not vice versa. This translates to a requirement 
for r such that rj • VT > pj 2 . 

We can also express the Fourier heat divergence in 
terms of reduced variables. 



j-V- 



. 1 -Idu. 



VT 



(23) 



Manipulation with Eq. 4 produces a form where the 
sign of u and directions of j and VT are irrelevant. 



Vq K 



3 



2ku 4 dT 



(24) 



Thus the sign of V • q K is determined by the sign of , 
which is valid for both p and n-type elements regardless 
of the sign of u. 

The Fourier heat divergence criterion is a convenient 
definition to distinguish these two regions of thermoelec- 
tric cooling in experimental data. The Peltier cooling 
region, defined by V • q K > 0, is found in the phase space 
where > 0. Likewise, the Thomson cooling region 

defined by V • q K < is the phase space where < 0. 
The constant relative current u(T) = const, separates 
the Thomson-type and from the Peltier-type solutions to 
the differential equation. In Figure 3, the CPM and u = s 
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cooler have opposite slopes, indicating these coolers ex- 
ist in separate regions of the cooling phase space. This 
result is consistent with our discussion above concerning 
the concavity of T(x) in Figure 5. 

For clarity, we suggest coolers which are predominately 
in the Thomson phase-space, V-q K < but may not have 
U = s be referred to as "Thomson coolers" . Similarly, 
"Peltier coolers" should refer to coolers operating in the 
usual V ■ q K > Fourier heat divergence phase-space 
where Joule heating dominates. 

This understanding of phase space for u = s and CPM 
coolers enables us to hypothesize that the performance 
advantages of u — s coolers extends to imperfect Thom- 
son coolers. We expect such coolers possess two primary 
advantages over traditional Peltier coolers. First, for a 
given material zT, performance (AT max and cf>) of the 
Thomson cooler is greater (Figure 2). The AT max solu- 
tion for the u = s cooler is compared to a Peltier cooler 
with the same material assumption for z in Fig. 4. Here, 
the maximum Seebeck coefficient is set by the band gap 
(Eg = 1 — 3 eV), per Eq. 19. The Thomson cooler pro- 
vides significantly higher AT max than the Peltier cooler 
with the same zT, nearly twice the AT max for E g = 3 eV . 
Second, in a Thomson cooler, the temperature minimum 
is not limited by zT explicitly like it is in a traditional 
Peltier coolers. 



6. DISCUSSION 

Efficiency improvements from staging and maintaining 
u = s also exists for thermoelectric generators, but the 
improvement is small (< 10% compared to CPM). This 
is because the u does not typically vary by more than a 
factor of two across the device. However, in a TEC the 
compatibility requirement is much more critical. When 
operating a TEC to maximum temperature difference, 
the temperature gradient varies from zero to very high 
values, which means u will have a much broader range 
(Figure 3) in a TEC than in a generator. Thus, unless 
compatibility is specifically considered, the poor compat- 
ibility will greatly reduce the performance of the thermo- 
electric cooler, and this results in the AT max limit well 
known for Peltier coolers. 

In real materials, changing material composition also 
changes zT so the effect of maximizing average zT is 
difficult to decouple from the effect of compatibility. As 
such, efforts which are focused on maximizing zT will 
generally fail to create a material with u = s and may 
only marginally increase AT max . Conversely, focusing on 
u = s without consideration of zT could rapidly lead to 
unrealistic materials requirements. 

In this new analysis we have focused on the compat- 
ibility criterion, u ~ s, with constant z (as opposed to 
zT [16]) to demonstrate the differences between a Thom- 
son and a Peltier cooler typically analyzed with the CPM 
model. Generally, achieving u — s in a material with fi- 
nite zT, is more important to achieve low temperature 



cooling than increasing zT. 

Minor improvements in thermoelectric cooling beyond 
increasing average zT by increasing the Thomson effect 
in a functionally graded material were predicted as early 
as 1960 [7]. Similarly Miiller et al. describe modest gains 
in cooling from functionally grading [8-10] where mate- 
rial properties are allowed to vary in a constrained way 
such that the average zT remains constant. Such addi- 
tional constraints can keep the analysis within the Peltier 
region, preventing a full optimization to a u = s solution. 

Bian et al, [11, 12] propose a thermoelectric cooler with 
significantly enhanced AT rnax using a rapidly changing 
Seebeck coefficient in at least one region. The method 
of Bian et al focuses on the redistribution of the Joule 
heat rather than a consideration of the Thomson heat 
or the effect of compatibility. Nevertheless, the region of 
rapidly changing Seebeck coefficient would also create a 
significant Thomson effect and likely place that segment 
of the cooler in the Thomson cooler phase space while 
other segments would function like a CPM Peltier cooler. 

In a traditional single-stage (or segmented) thermo- 
electric device, the current flow and the heat flow are 
collinear and flow through the same length and cross- 
sectional area of thermoelement. This leads to the com- 
patibility requirement between the optimal current den- 
sity and optimal heat flux to achieve optimal efficiency. 
In a multi-stage (cascaded) device the thermal and elec- 
trical circuits become independent and so the compati- 
bility requirement is avoided between stages. 

A transverse Peltier cooler also decouples the current 
and heat flow by having them transport in perpendicular 
directions. Gudkin showed theoretically that the trans- 
verse thermoelectric cooler, could function as an infinite 
cascade by the appropriate geometrical shaping of the 
thermoelement [27]. The different directions of heat and 
current flow enable, in principle, an adjustment of geom- 
etry to keep both heat and electric current flow indepen- 
dently optimized. Cooling of 23 K using a rectangular 
block was increased to 35 K using a trapezoidal cross- 
section [28]. 



7. CONCLUSION 

Here, we compare self-compatible coolers with CPM 
and commercial Bi2Te3 thermoelectric coolers. Signif- 
icant improvements in cooling efficiency and maximum 
cooling are achieved for equivalent z when the cooler is 
self-compatible. Such improvement is most pronounced 
when the goal is to achieve maximum temperature dif- 
ference, rather than high coefficient of performance at 
small temperature difference. Optimum material profiles 
are derived for self-compatible Thomson coolers and re- 
alistic material constraints are used to bound the per- 
formance. Self-compatible coolers are found to operate 
in a fundamentally distinct phase space from traditional 
Peltier coolers. The Fourier heat divergence of Thomson 
coolers is dominated by the Thomson effect, while this 
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divergence in Peltier coolers is of the opposite sign, in- 
dicating the Joule heating is the dominant effect. This 
analysis opens a new strategy for solid state cooling and 
creates new challenges for material optimization. 
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By conservation of energy (or recursive relation Eq 25), 
the heat exiting the n pumps, Q n is the heat pumped by 
the first pump plus the total power, Eq 29. 



Qn — Qo + P 



(31) 



Combining equations 30, 31, and 28, it is straightfor- 
ward to show 



1 n 1 



(32) 



8. APPENDIX 

The metric for summing the efficiency and coefficient 
of performance of these thermodynamic processes is not 
a simple summation because energy is continuously be- 
ing supplied or removed so that neither the heat nor en- 
ergy flow is a constant. The derivation, attributed to 
Zener[] 7], of the coefficient of the performance Eq. 2 and 
efficiency summation metric for a continuous system in 
one-dimension given here is based on Zener and similar 
derivations given in [7] [29] [30] [19]. 

Consider n heat pumps (or heat engines) connected in 
series such that the heat entering the ith pump, Qj, is the 
same as the heat exiting the i — 1 pump, namely Qi-\. 
Then by conservation of energy 



Pi 



(25) 



where Pi is the power entering the pump i. The coeffi- 
cient of performance <p n of the pump n is defined by 



1 = Pi 

4>i Qi-1 

This gives the recursive relation 

Qi = <9i-i(i + x) 

which can be solved for Q n as 

n 

Qn = Qo~[[(l+-) 



(26) 



(27) 



(28) 



The total power added to the system P from the n 
pumps is 



(29) 



i=i 



The coefficient of performance for the n pumps, <p, with 
heat Q c = Qo entering from the cold side is defined by 



1 _ P 

4> Qo 



(30) 



This can be transformed into a summation by use of a 
natural logarithm 



1 ™ 1 

>(i + T) = £Ml + -r) 



(33) 



i=l 



While 1 /<pi should become very small with small 5Ti 
Ti — Ti-i the reduced coefficient of performance (p 
should remain finite 



)_SJ\ 



(34) 



Assuming a monotonic temperature distribution (for a 
simple TEC Tj > Xi_i for all i) the sum in Eq. 33 can be 
converted to an integral in the limit that n — > oo, where 
STi = (Th — T c )/n and ln(l + x) — > x for small x 



ln{l + \) 



1 



T c T(pr{T) 



dT 



(35) 



If the temperature distribution is not monotonic but can 
be divided up into monotonic segments, each of these 
monotonic segments can be individually transformed into 
integrals. 

For a generator (as opposed to a cooler or heat pump) 
power is extracted rather than added at each segment 
in the series. Then the sign of the power in equation 
25 is negative for a generator with the efficiency given 
by i] = —P/Qa = —l/4> and reduced efficiency r\ r — 
— l/(p r . Thus the above method can be used to derive 
the analogous equation for generator efficiency: 



ln(l — i]) 



T ' 



dT 



(36) 
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